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1. INTRODUCTION 

The concept of networked control systems (NCSs) have been mentioned by many scientists in their 
works as early as the end of the twentieth century. Because the low maintenance cost and convenient installa- 
tion, NCSs have aroused widely concerned. Due to the introduction of shared network, new constraints occur 
when the plant outputs and control inputs are transmitted through communication networks: quantization errors 
in the signals transmitted through the network, packet dropouts, network-induced delay, access constraints and 
power consumption mainly in wireless networked control systems [1]-[5], thus increasing the complexity of 
the system. These factors will affect the reliability of the system, and cause the system performance decline. 
In order to improve the reliability and security of the NCSs, fault diagnosis has been widely used in engineer- 
ing systems such as aero engines, dynamic vehicle systems, chemical processes, and power systems [6]-[8]. 
Fault detection (FD) has received widespread attention as one of the most considerable parts of fault diagnosis. 
Failure is the phenomenon that the state of the system deviates from the normal working range due to satura- 
tion, stuck or degradation of actuators, sensors and other components, which has a negative influence on the 
system performance. As a result, it is very significant to detect the system faults as soon as possible to ensure 
the safety of systems. There are many researching results on FD for kinds of systems with various methods 
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[9]-[12]. Pan and Yang [11], Hə FD filter is proposed for NCS with large transfer delays. An adaptive Kalman 
filter for NCSs is proposed in [12] to minimize the effects of delay on the residual signal. The fault detection 
system for NCSs with packet losses was been designed by modelling the NCSs as a Markov jumping linear 
system [13], [14]. For research work of FD problem of NCSs, construct appropriate filters and state observers 
as residual generators to generate residual signals [15]-[17]. The problem of FD of a NCSs under communica- 
tion constraints limited is considered [18]. The practical guideline for selecting the right protocol in industrial 
networked control systems (NCSs) is provided in [12], [19]. The evaluated the performance of MPOLSR pro- 
tocol and MDART protocol using NS-2 based on the success delivery rate and packet loss bas been studied in 
[20]. Inspired by the above discussions, the main goal of this paper is to devoted to FD of NCS subject to both 
random communication delays, stochastic packet losses, limited communication and access constraints. The 
main contributions of this work can be highlighted as follows: 


- The fault detection problem is extended for a class of networked control systems (NCSs) with random packet 
losses, time-varying delays and limited communication to reflect more realistic environment, 


- Adaptive Sliding mode observer approach is utilized to deal with the fault detection, 


- Residuals generator is designed, reducing the false alarm rate. Then, the residual signals are evaluated and 
compared with a threshold to detect the faults occurrences. 


- Application to a Hydro-turbine governing system [21], [22] shows that the proposed method achieves better 
fault detection. 


- True-Time toolbox is used to reflect a more realistic numerical network communication and validity of the 
proposed design method. 


- The control performance of the proposed method is evaluated for several industrial protocols: EtherNet/IP 
protocol, PROFINET IRT protocol and EtherCAT protocol of the standard IEEE 802.3, the right protocol for 
a NCS is provided . 


The remainder of this paper is organized as follows: section 2 introduces the problem statement and 
preliminaries, our step of adaptive observer synthesis is given in section 3. The simulation results based on 
True-Time toolbox and Matlab/Simulink will be given in section 4 to verify the efficiency of proposed method. 
Finally, a conclusion is provided, including some perspectives of this work. 


2. PROBLEM FORMULATION AND PRELIMINARIES 
In this paper, the discrete linear system with output delay is structured as Figure 1; the state-space 
model of the linear plant dynamics (1). 


a(k +1) = Ax(k) + 2,2(k — Tk) + Bu(k) + Td(k) + FY(k) 
y(k) = Ca(k) (1) 


Where x(k) € R” denotes the state vector, x(k — Tk) € IR” denotes the state delay vector , u(k) € R™ denotes 
the control input vector Y (k) € R4 is the fault vector, y(k) € R? denotes the measured output vector and 
d(k) € R” the noise vector , A , 2,, , B , C and F are matrices of appropriate dimensions. 


Assumption 1 [13]: It is supposed that random packet losses exists in output channel. It is modelled in 
the system as Bernoulli process. We define y(k) the output of the system (with internal noise) and y(k) the 
data used with a probability 8(S, {8, = 1} = 3). If the data is not available, we will use the preceding data 
y(k — 1) with probability 1 — 8(S, {8, = 0} = 1 — 8). 

The following equations describe this phenomenon. 


y(k) = By(k) + (1 — B)y(k — 1) (2) 

Where 6; € {0,1} obeys the Bernoulli distribution. 
Assumption 2 [23]: In this paper we will consider that the band-width of the communication network 
connecting the sensors and the fault detection module which generates a residue is limited capacity, w, sensors 


among p can reach these channels to communicate with the residues generator while the others remain on 
standby. Similarly, wọ from p actuators receive their command from controller at each sampling period. 
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l < wo <m 
l < w; <p (3) 


Taking the phenomenon of sensor and actuator saturation, the functions of saturation ¢(-) : R” > R” 
and o(-) : R” — R” are defined as (4). 


c (k) = [cP (k), oF (k), .sh(h)] 
o(k) = [01 (k), 03 (k), , 0n,(k) 


(4) 


\" 


If we denote by y(k) denotes measurement vector available to controller and fault detection module 
and u(k) denote the control signals generated by controller at discrete time k. Based on the above communi- 
cation sequence, We get the relations in (5). 


M.(k).y(k), M; = diag (si(k)) (5) 
Mo(k).u(k), My = diag (0:(k)) 


y(k) 
u(k) 


From (1), (2), (3), (4) and (5), the dynamics of the networked control system via a shared communi- 
cation medium can be described as (6). 


x(k +1) = Ar(k) +2,.0( k= Tk) + BM,(t)julk)+Trd(k)+ FYT(k), 
y(k) = By(k) + (1 — B)y(k — 1) (6) 
ylk) = M.(t)Cx(k) 


Lemma 1 [24] (Schur complement). Given constant matrices of appropriate dimensions %11, B12 and Baz € 
R”"*”, where B1 = Bi. Boo = Bs, then Bə > 0, 811 — $1.85, BT, > 0 if and only if 


Lemma 2 [25] Given matrices of appropriate dimensions =11 = =e , 212 and =99. § a function which satisfies 
¥ = 8 < I, where is an identity matrix. Then the inequality =11 + 2126222 + ean < 0. Is not true 
that, if and only if there exists a scalar 3 such as the inequality is checked: =), + 3212 =e + 37l =e < 0. 
or equivalently 


Sar Je 59 
* —3F 10 0 <0 
x x —3=12 


where the symbols (*) denote the symmetric terms. 

The remote FD dynamic behaviour of NCSs is illustrated in Figure 1. Notably ,y(k) denotes the ac- 
tual output signal and y(k) denotes the output signal used by controller,u(k) is the control signal produced 
by controller and u(k) is the actual control input. In this structure the fault information is not affected by the 
communication delay between the sensor and the node of FD unit and communication delay between controller 
and the node of FD unit. 


Performance evaluation of industrial ethernet protocols for real-time fault detection ... (Samba Aimé Hervé ) 


264 g ISSN: 2089-4856 












y(k — Tk) ii(k) 


Communication 
network 


By = 0 y(k) 





y(k) 





Fault 
Detection Unit 


Figure 1. The proposed block diagram of remote FD in networked control systems with time delays, packet 
losses and access constraints 


3. RESIDUAL DYNAMIC SYSTEMS 

In this section, we aim to design an fault observer based on discrete-time adaptive sliding mode for 
considered networked system with packet losses, communication delays and access constraints. The structure 
of the adaptive fault observer proposed is given by (7). 


t(k+1) = Afk) + 2, 2(k — Tk) + BM,(k)u(k) + Lly(k) — ilk) 
T (k +1) = AT (k) + Ao (y(k) — G(A)) (7) 


glk) = (1 — B)M.(k)C2(k) 


Where ĉ € R” and (k) € R? denote the state estimation for the fault observer and estimation of the 
measurement output respectively; L € IR”*? is the the observer gain matrix and Y (k) is the estimation of the 
fault vector. This adaptive sliding mode observer allows to generate a residue which will be analyzed to detect 
the faults. It is significant to note that in the context considered in this paper, A1, Az are the fault detection 
parameters. 

Let €,(k) = x(k) — &(k) and e,(k) = y(k) — ĝ(k). 

According to (1) and (7), we get the estimaror dynamics by (8). 
Es(k +1) =(A— (1 — B)M,LC) e(k) + (1 — 8) LCax(k) + Td(k) (8) 
+ FY (k) == A- Erlk = Ty) 


Ey(k) = (1 — 8)M,(k)Cex(k) (9) 
where Y (k) = Y(k) — T (k) 


For the fault observer, the logic of fault detection considered in this work is given by (10). 


3(k)|| > Ten Alarm, fault is detected (10) 
ISk) < Sen No alarm, fault is no detected 


Where J(k) is the residual evaluation function of the residual generator and Jt, is the threshold are 
selected as (11). 


Sth = sup KE {3(k) } (11) 
d(k)€I2T(k)=0 


where n is the length of the evaluation window and r(k)) = ¢,(k) denotes the residual signal of the system. 
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4. STABILITY ANALYSIS 


Theorem 1. Consider the system (7) with the given scalar y > 0, if there exist matrices P; > 0, Aj, 
A> > 0 and L, such that the following matrix inequality holds. 


211 yà X 


H= * —y A2 0 <0 (12) 
‘i t = 19 
where 
= | -Pi + 1 0 
Eu =| 0 aoe 
—P 0 
X22 = | 0 i E | 
Sn Ai 
Jas A z 
j | X21 299 | 
~ = T 
Ži = (Ao(1 — B)M.C) 
Sy = FTP, 
d/99 = Ad 


Then the residual dynamic system (7) is asymptotically stable. 
Proof 1. Choose a Lyapunov functional candidate as (13). 


V(k) = El (k)Piex(k) + YT (k) Y (k) (13) 
Then , it can be obtained as (14). 


E {V(k +1) — V(k)} =E{e2(k)? (A — (1 — B)M-LC)* Pı (A — (1 — 8) MLC) es (k) 

A — (1 — )M,LC)" PiLFÝ(k) + ÝT (k) FTP FÝ(k) 

EYT (kIT (k + Des(k =r) AL PiN erlk — Te) 

Denk) (A158) MLC KAn L erlk = Th) (14) 
+ 2e7 (k)(A = (1— BM LOV (A—(1— B)M.LC)' Pi 

+ 2e(k — Th) A,,PiPd(k) + [Pd(k)]* Py [Pd(k)] 

— €,(k)? Pyex(k) — 7 (k)Y(k) 


where 


© 
Hc 
> 
+ 
cee 
KH 
| 


-A T (k) — Aves (k) + &(k) 
Plk) = A, T(K) + Te +1) (15) 


Taking ||d(k)|| < d, where d is know positive constants, it is derived that 


[Pa(k)|* Py [Pa(e)] < |Pa(e)|? Pall < [IP a)” Ps I (16) 
[A — (1 — B)M.LC] *P, [C'd(k)] < || [A — (1 — 8) MLC] *P,|| || a(k)|| (17) 
2e2(k — Tk)” AL Pr d(k)] < 2 |len(k — Te)" A || Pal] [IE] d] (18) 
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According to (15), (16), (17) and (18). It can be further obtained (19). 


E{V(k +1) — V(k)} k) {(A— (pM Pi (A — (1 — pM) — Py} 
isu hee Qex(h AAS alas + TT (k) 
x (FTP T(k) —DY(k) + 2||(A -— (1 — B) MLC)" P: || 


x |lex(k)” |IIVI dl] + 2 |le(& — 7) a . || IP 
x [Tl] d] + e2(k)" (A- (1-8) MLC)" Py 


x (A — (1 — B)M.LC)e,(k) + 2\762(k)? Pree (k) 


+MTOT + &(k) (MTO + 8(k)) Pill + [IPI] 


where, 


According to the Lemma 2 and (19), it easy to obtain (20). 
E{V(k +1) — V(k)} < MHM + 2M’ OT 6(k) + &(k)* (k) 
For the matrix H, it can be obtained by Lemma 1: 


Ait FA X 


H = x —y 22 0 <0 
= i =Y 
where 
= — P; ae yI 0 
A | 0 -I +I | 
—P 0 
2192 = | “I | 
Sur Nye | 
A= & = ; 
a | dig, 299 


~ = T 
i= (A —(1- B)M,LC) Pis 
~ = T 

2 (A2(1—- B)M.C), 

Xin) = F' Py, 

ding = At, 

The proof is completed. 


(19) 


(20) 


(21) 


Theorem 2. Consider the dynamic system (8). If there exist matrices L € R”*™ and 7? < 1 


satisfying the condition. 


orn Oyo 
| * No | < 0 
Where 
Qi = =n" G = )C) * ((1— B)M,(k)C), 
Qiz =2(A- BM IO). 
se = aad ((1 — B)M,(K)C), 


2 
then system motion gets into the sliding surface in finite time. 


(22) 
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Proof 2. According to [22], the system motion gets into the sliding surface in finite time, then s(k) 


satisfies. 
ls(k + 1)I| < |Is(*)| (23) 
the condition (23) can then be reformulated in (24). 
s(k + 1)" s(k + 1) — n?s(k)?9(k) = ef (k + Ley(k + 1) — ne? (k)ey(k) (24) 
According to (8), then we have 
s(k +1)’ s(k +1) — n°s(k)’ s(k) =e; (k + 1)((1- - <(k)C)* ((1 — B)M,(k)C) 
X E(k +1) — n*er (k)((1 — B)Mo(k)C)” (25) 
x ((1 — 8)M. P )C)Ez(k) <0 


where 77 < 1, then substituting (8) into (26), we can obtain (26). 
s(k + 1)" s(k + 1) —°s(k)" s(k) =e; (&)[(A — (1 — b)M LC)" ((1 — 8) M.(k)C)” 
x ((1 — B)M.(k)C)(A — (1 — 8) MLC) 
-= n? ((1 — B)M.(k)C)" (1 — 8) M.(k)C) 
xe (kite k= Th) AF ((1— B)M.(k)C)* 
x ((1 — B)M,(k)C)eu(k — Te) + [Pa(k)]" (26) 
x ((1 — B)M.(k)C)" (1 — B)M.(k)C) 


LTdlky t21 (hk (A= C= MLC 
x ((1 — 8)M,(k)C)" ((1 — 8)M.(k)C) An, 
x €,(k — Tk) + 2e} (k) (A — (1 — B)M.LC)~ 
x ((1 — 8)M,(k)C)" ((1 — 8)M,(k)C)P'd(k) 


+ 2e (k) (A — (1 — B)M; LO) AT 
x ((1 — B)M;(k)C)" (1 — B)M,(k)C)Id(k) 


Simolarly, (24) can be obtained (27). 


skti s(k + Dnr s(k) ce (k)[(A — MLC (1 — Mk) 
— n’ ((1 — B)M,(k)C)* ((1 — B)M,(k)C 

A—(1—B)M,LC)((1 — 8)M,(k)C)len(k)) 

2e7 (k)(A — (1 — B) MLC)" (1 — B)M,(k)C)* 

((1 — B)M.(k)C)L7,€n(k — Th) — 2e} (k) 

(A T (1 = B)M,LC)* ((1 p B)M,(k)C)* 

((1 — B)M.(k)C)Ao(k) — 2en(k — Tk)” (27) 

((1 — 8)M.(k)C)* ((1 — B)M,(k)C) A(R) 


x 
+ 
x 
x 
x 
x 


One can further get (28). 
s(k +1)" s(k + 1) — °s(k)* s(k) <ez (k)[2(A — (1 — 8) MLC) ((1 — 8)M.(k)C)* 
x ((1 — 8)M,(k)C)(A — (1 — 8) MLC) (28) 
— n’ ((1 — B)M.(k)C)* ((1 — B)M,(k)C)| 
xX Ex(k) < 0 
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According to the Lemma 2, we obtain (29). 


s(k +1)’ s(k +1) -—n’s(k)" s(k) < ef (k)Zes(k) < 0 (29) 


This completes the proof. 


5. SIMULATION RESULTS 

In this section, we will propose a numerical example of simulation to illustrate the effectiveness of 
the methods presented in this work. Let us consider the model of networked control hydroelectric power plant 
[20]. The flowchart of the winnowing device control and communication network is represented on Figure 2. 









y(k — Tk) i(k) 






Communication 
network 


B, -0 | I) 


Figure 2. Schematic diagram of networked control hydroelectric power plant 


This system is been used in [21], [25], for the design A Networked iterative learning fault Diagnosis 
algorithm for systems with sensor random packet losses, time-varying delays, limited communication and 
actuator failure. 

The state representation of dynamic model is described as (30). 


1.1840 —0.4046 0 1 
x(k+1)= | 05000 o0 O|a(k) +] 0 | u(k) 

0 0.5000 0 0 (20) 
y(k) = | 0.2943 0.3382 0.0001 | (k) 


0.034 0 —0.01 
Define A,, = | 0.031 0.03 0 
0.04 0.05 —0.01 
According to the time scale, we define the time-varying communication delays as 7;(k)(z = 0,1,2). 
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Let the fault signal Y (k) be given as (31). 


; .2sin(k) 122 < k < 142 
rw=f 0.5 + 0.2sin(k) <k < 31) 


0, others 
AndF=[1 11 


The communication constraints is fixed to one channel (w; = w, = 1) so the following 3-periodic 
sequence can be proposed. 


SoS oo 
GS SS 
roo Tee 
Soja ee 


The random variable 8 satisfies the Bernoulli distribution, let 8 = 0.65. 
Applying Theorem 1 and 2, we can obtain according to the Matlab LMI toolbox, the desired P , the 
observer gain L, the scalar 7 and parameters A, and A, as follows: 


0.231 —O0.005 0.002 


Py = 0 0.0012 0.0021 |, 
0 0.0321 0.0063 
—0.0352 
L= | 0.000011 | , A, = 0.03621 , A2 = 0.0073 and 7 = 1.0173e7°. 
—0.0142 


We obtain the Simulink/ true time models of networked control hydroelectric power plant Figure 3. 
The true-time network block simulates the access to the medium and allows the transmission and the reception 
of data through the network. Table 1 shows the simulation parameter for wireless network block. 


Table 1. Simulation parameter for network block 


Parameter Values 
Network type 802.15 (LAN) 
Data rate 200 Mps 
Minimum frame size 544 bits 
Transmit power 200 dbm 
Receiver signal threshold -48 dbm 
Path loss exponent 33.5 


In order to show the effectiveness of the approach proposed, the noise signal d(k) is a white noise 
Gaussian of an amplitude of 0.025 (sample time T; = 1s). The fault signal Y (k) occurs between the moments 
122°” at 157°” steps. 

The results comprise a packets loss at the moments: steps 11 at steps 19 steps, steps 41 at steps 52, 
steps 63 at steps 68, steps 84 at steps 86, steps 89, steps 101 at steps 108, steps 128, steps 133 at steps 139 and 
steps 194 at steps 196 Figure 4. 

Figure 5 show the evolution in real time of the actual speed of hydroturbine and the estimated speed by 
adaptive sliding mode observer with EtherNet/IP protocol, PROFINET IRT protocol and EtherCAT protocol. 

The generated residue is illustrated by the Figure 6, which shows that the residual signal converges to 
zero without faults Figure 7 , then , changes rapidly when the faults occurred this for the three communications 
protocols. Figure 8 illustrate the residual evaluation function J(k). 

Figure 9 shows that the same faults is detected at steps 125 with switched ethernet protocol. Figure 
10 shows the occurrence of fault can detect the fault at steps 124 with EtherCAT protocol. 
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Figure 3. Simulink/true-time model of remote FD in networked control hydroelectric power plant 
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Figure 5. Actual and estimated speed of hydroturbine 


Int J Rob & Autom, Vol. 10, No. 3, September 2021 : 261 — 274 


Time in sample (k) i 


Oe 


200 


Figure 4. The distribution of packet losses ” 1” means packet received, ”0” means packet lost 
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Figure 7. Residual signal r(k) without fault 
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Figure 8. The residual evaluation function J(k) 
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Figure 9. The residual evaluation function §(/) with EtherCAT protocol 
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According to (11), we can obtain the threshold of fault detection for EtherNet/IP protocol, PROFINET 
IRT protocol and EtherCAT protocol . The values of J:n and time steps that fault can be detected are reported 
in Table 2. We can notice that when fault occurs at steps 122, the occurrence of fault is detected at the steps 
126 with EtherNet/IP protocol which corresponds to a delay of detection of steps 4 Figure 11. 

According to Table 2, thecarriersensemultipleaccess /collisiondetection EtherCAT protocol will 
be adapted better to improve the performances of faults detection module in networked control hydroelectric 
power plant. 
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Figure 10. The residual evaluation function 3(/) with PROFINET IRT protocol 
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Figure 11. The residual evaluation function J(&) with EtherNet/IP protocol 


Table 2. Threshold +n and time steps for different protocols 


Protocols Threshold (Jn) Time detection 
EtherNet/IP 0.04919 steps 126 
PROFINET IRT 0.0485 steps 125 
EtherC AT 0.0524 steps 124 


6. CONCLUSION 


In this paper, an adaptive sliding mode observer-based fault detection for a class of NCSs subject to 
random communication delays, stochastic packet losses, access constraints and bounded disturbances has been 
established. A Lyapunov stability theory based on stochastic function is presented to analysis the asymptotical 
stability of the fault detection systems. The dynamic hydroelectric power plant model is used to verify the 
effectiveness of the proposed method based on TrueTime and Matlab/ Simulink, corroborated our predictions 
that a ethernet for control automation technology (EtherCAT) protocol would be more appropriate to reduce 
the false alarm rate and increasing the efficiency of the remote control hydroelectric power plant. Furthermore, 
the proposed results can be extended to the FD of NCSs with random communication delays, stochastic packet 
losses, access constraints quantization effects and possibly time-varying uncertainties. 
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